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SUMMARY 

The application of the well-known hasic principle of 
mechanics, the principle of Jourdain, to pro'blemB of the 
theory of the boundary layer ^le^ds to an equation from 
which the equations of Von Karman, Leibenson, and Golubev 
are derived as special cases< The given equation may he 
employed in other integral methods. The present paper 
deals with the method of the variation of the thickness 
of the "boundary layer. A number of new approximate for- 
mulas valuable in aerodynamic calculations for the fric- 
tion distribution are derived from this procedure. The 
method dias been applied only to laminar boundary layers, 
but it seems probable that it may be generalized to in- 
clude turbulent layers as well. 



The first step in the field of application of inte- 
gral method^ to the theory of the boundary layer was made 
by Von KaVman in 1981. The fundamental idea of the 
Von Karman method (reference l) consists in replacing 
the true velocity distribution in the boundary layer by 
certain approximate velocity distributions satisfying the 
differential equations of the boundary la^er only on the 
average. Tor the averaging method Von Karman applied the 
usual integral mean over the layer, a procedure which led 
to a certain integral condition which is no other than 
the theorem on the change of momentum of a fluid in a 
volume element of the boundary layer. In the application 
of this method fundamental significance is attached, of 
course, to the choice of the system of functions that 
replaced the true velocity distribution, a circumstance 
analogous to the corresponding methods of the theory of 
el^asticlty (inethods of Ritz , ftalerkin, and others) and 
presenting, in general", great difficulties in using the t. 
so-called "direct" methods. 



*Inst. of Mechanics of the USSR Academy of Sciences, 
Applied iiathemfttics and Mechanics , vol, V, nOf 3,1941^ 
pp. 453-470. 
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In attempting to enjploy polynomi^al ^approximations 
for the velocity di st r i"but ion , Von Karman found it was 
necesBary to diapens© with the theory of the "boundary 
layer of infinite thickness (asymptotic approach of the 
velocity in tfee layer to the velocity of the outer flow) 
and introduce the concept of the layer of finite t'fclck— 
ness. If, instead of pclynomials, transcendental inunc- 
tions are applied as the approximating fui^ctlons, the 
usual assumption croncerning the asymptotic boundary 
layer can he used. 

The method of Von Kayman came into general use and 
was successfully applied in practice. At the same time, 
hovfevQTf a numter of investigators have pointed out cer- 
tain defects: namely » the restricted field of applica- 
tion and the ahsence of any real basis (reference 2). 
The existence of a single integral condition permits 
making use Of a family of functions approximating the 
velocity profile in the layer with the functions having 
only one parameters namely, the thickness of the "bound- 
ary layer* The approximation i;B obtain-ed by making use 
of the boundary conditions at the wall and at the edge of 
the layer. The number of boundary conditions, both those 
given in advance and those derived from the equations of 
motion, is very limited; and, moreover, among these condi- 
tions the value of the first derivative of the velocity 
along the normal to the surface of the body does not 
enter. This derivative is precisely the fundamental 
unknown proportional to the frictional stress at a given 
point of the body. The choice of any particular boundary 
conditions leads to "ejxternal" solutions more accurately 
describing the phenomena near the outside edge of the 
layer or to "internal" solutions approxtmAting the flow 
conditions near the walls, (See reference "S, ) The well- 
known solution of Pohlhausen (reference 4) is an exter- 
nal = solution. In setting up the approximating poly- 
nomial of the fourth degree three conditions were employed 
(the values of the function and of the first and second 
derivatives) at the outer edge and only two at the wall, 
the first derivative being of necessity omitted from the 
other two conditions^ only the values of the fu^iotion 
and of its second derivative being used. The Pohlhausen 
solution gives an exaggerated friction and a retarded 
boundary layer separation-. At the jsuthor^s suggest ion, 
A. m. Alexandrov (reference 5) applied a second system of 
conditionsf leading to an internal solution, a reduced 
friction, and a some^fhat earlier separation. By making 
use of a polynomial of the sixth degree a cert^-ln added 
accuracy could have been obtained. 
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Further progroBS in the development of integral 
methods w?i8 made by setting Tip equations which were a' 
gc^ne:^ali«ation of the momentuffl equation given hy Von 
Earman. L. S, Leihenson (reference 6) proposed as a 
second integral condition the equation of the change of 
the kinetic energy of an element of the boundary layer, 
the condition being obtained by integration of both sides 
of the Prandtl equation, first multiplied by the forward 
velocity along the normal to the wall from the surface 
of the body to the edge of the layer, V. V. G-olubev 
(reference 7) pointed out the infinite possibility of 
other integral conditions obtained by the sajn© device 
of mulliip lying the two sidea of the Prandtl equation by 
successive integral powers of the forward velocity, The 
same method for the paarticular case of a flat plate was 
followed by 7, Sutton (reference 8)",. Kone of these in- 
vestigators made any attempt to indicate methods for mak~ 
ing -use of the derived new integral conditions. Only 

Sutton gpve an example of the computation of the fric- 
tion of a plate,.making use of a set of two equations. 
It is not clear, however, whether any of the proposed 
conditions could be nsed individually or whether it was 
necessary to use two of them combined. 

In the present paper use is made of a single general 
principle of mechanics which, as far as is known, has not 
yet been applied in hydrodyanraic s and, in analytical me— 
ch»nics bears the name of the Jourdain principle (refer- 
ence 9). This principle, which is intermediate between 
the well— known principles of D'Alembert and G-auss, is very 
convenient for problems of hydrodynamics in the applica- 
tions of Euler and, in part icular ,; for obtaining "direct" 
methods of the solution of boundary layer problems, 

Ih© assigning of a particular form for the velocity 
variation in the expression of the Jourdain principle leads 
to the equation of Von KfTrman, the system of equations 
of Leibenson— frolubev , and other methods which may be use- 
ful in applications. 

1. The principle of Jourdain is obtained from the 
principle of D'Alembert by differentiation witja respect 
to time. It is convenient to take the general equation 
of- mechanics in th© forjaj 



n 




« 



(1.1) 
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whero , vri , and Jj^ ftre, respectively, mass, accel- 
erRtiony and giveaforce for a. certain 1th point and . 6rj^ 
Is n virtual displacement. By differentiating with re- 
spect to time, there is obtainedJ 

n 



1 



n 

dSri 



(miWi - Pj) + y -^-CmiWi I'i)firi 0 (l,3) 

dt it ^ 

i=i i=i 



By the definition of virtual displaceraentes 



dt ^ 1 1 



Traene is the velocity of the ith point, 

let the actual motion of the system at a given 
instant be compared with a neighTsoring motion differing 
from^ the true one in the velocities but not In the posi- 
tions - that is, let 6ri « 0. Equation (1.2) then 
assumes the f-orm: 

n 

^ (m^Wi - I'i)Bvi =s 0 (1.3) 

This is the general equation of the principle of Jourdain. 
It is to "be enphasiaed that are given forces while 

the reaotione of the idfeal connections % satisfying 
the condition 



n 



i«x 



1% 6rj_ » 0 



do not enter equation (1,3), since "by the fpregoing ex- 
planation 



^i = 0 

The fricttonal forces, as always, 1i^sl«ag to the class of 
given forces. 



In applying the Jourdain principle to the motion of 
a viscous Incompressible fluid the equation of the prln- 



I 
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ciple may "be written in the following most general form: 



where t is an artiitrary volnine of the fluid "bounded by 
the surface crj Sv is not suhject to any reatrictions 
except the condition of the existence of the integral on 
the left side of the eq^uation. 

It may he noted that, on the part of the surface cr 
corresponding to a solid wall, the variation Sv should 
not, in advance, he suhject to the condition (8v)j^ « 0 

since normal pressures enter the expression for the prin- 
ciple. In the same way, from the condition of the adher- 
ence of the fluid to the wall it does not follow that 5v 
must he set equal to zero because in the equation (1,4) 
the fractional forces are taken into account. 

If it is assumed that all the forces, hoth given and 
reactions, constitute an independent system, the fluid 
may he regarded as a free system of points not subject to 
any ri^lations. The condition of the conservation of mass 
for a constant density div v = 0 likewise imposes no 
restrictions on the variations of the velocity. Various 
particular principles and theorems of mechanics will h© 
obtained by restricting the variation of the velocity. 
Thus, for example, the momentum theorem is ob^tained by 
setting 6v = where € is an infinitely small vector 
independent of the coordinates; the theorem on the change 
in the kinetic energy is obtained by setting 6v = vc 
where € is an infinitely small scalar independent of 
the coordinates; and so forth* 

Considering only the case of steady motion and volu- 
metric forces having a potential Tt the principle is 
rewritten in the form 





pw - pF + grad p - {xV ^ 



v)Sv dT = 0 



(1.4) 




X V + grad B + vrot rot v) 6v $T = O 



(1.5) 



where 




whence by simple vector transformations there is obtained 
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(rot ^ X. v)6v dT * dfv(B5v)dT - /b div 8v dT 

+ V j iiv(rot V X 5v):cLT + v y rot v X 5 rot v dT s= 0 (1.6) 

U the second and fourth Integrals, in passing from volnae 
to surface integraXa* there is obtained 

[rot V X v)6v dT+ BSVjj da " ^ • 
+ V (rot V X 8v)a dpr + - v S rot v X rot v dT a 0 (1.7) 

The velocity variation 5v now is stil>^ee.ted to the two 

conditionsi 

1. Ihe neighboring motion satisfies the condition of 

incompr^lasibility div 6v = 0, 

2. On the boundary surface a the velocity variation 

6v is equal to zero. 

Unjder these assumptions the fundamental equation (1,7) 
assumes the form:' 

/ (rot vXv) 5v dT4- |5 /rot v x rot v d T « 0 . (1.8) 

This is the variational form of the principle of Jourdaln 
in hydrodynamics for very general restrictions imposed dn 
the velocity vatiation. 

For particular cases of motioii the prineHp^e takes a 
s t i 1 1 s'imp 1 er ' f orm, 0 ons i dOr •, f or eacampl e , the gt eady 
motion of a viscoUfp, inconvpressible fluid at very small 
Reynolds numbers. The inertia terms in the hydrodynamic 
equations may be neglected. In equation (l<.8) the first 
integral drops out and this leaves 
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<1.9) 



which la no other than the known principle of the mininmm 
dissipated energy first enunciated "by Helmholte and 
Eayleigh, Equation (1,9) maintains the same simple form 
also for any Eeynolds numhers if the motion occurs ^srith 
vortices parallel to the velocity ("free" vorticas of the 
wing). 

If equation (1,8) Is rewritten in the equivalent form 

rot V X rot V dT* O (l.lO) 




It is noted that, 'by suljjecting th? velocity variation to 
still another restriction, 



V X 6v s= 0 



(1»11) 



that is, if it is required that the varied motioa have the 
same direction of the velpoitle? at a given instant as the 
tru0 velocities, there is again oTstained the principle in 
the simple form (1,9), The general form of the velocity 
variation now will he 



6v ve<x, y, «) (1. 12) 

where the infinitely small function e(x, y, z) is arhi- 
trkry and is suh^ect only to the restriction . 

div 5v = £ div V + V X grade = 0 

which, for an Inoompresslhle fluid, is transformed into 
the condition,* / . " 

V X grad € - 0 (l. 13) 

This condition has a simple meaning: namely, the stream- 
lines should lie on the level surfaces of the functioR..€, 
In trro-dimenslonal motion with a*ial symmetry this simply 
means that c should he an arhitrary function of the 
stream function \|/ .— that Is, 
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6v = V e (vl/ ) (l. 14) 

The condition that the velocity variation TDOcome zero 
at the "boundary a reduces to the condition that the 
closed surface a should "be a surface of flow. 

In application, it is found simpler to make use of 
the principle in its general form (l,4), since the impos- 
ing of the a"bove conditions on the velocity variation 
involves great difficulties. Application of the principle 
in the boundary layer theory will now "be considered. 

Sju-TI|g_£age^of \a. tgOT-dimensional^ stead y l aminar 
boundary .layer , of _f initethickneis.— The usual considera- 
tions of the Prandtl theory'wlth yegard to the relative 
smallness, for largQ Beynolds num'bers, of the transverse 
lengths and velocities lead to the principle of Jourdain 
in the form 

■/ \ by p dx Sy2 / 

where a is an ar"bitrary region of the "boundary layer, 

Chooaing for the region a the part of the layer bounded 

by the straight lines x = , x = Xs , the wall of the 

body and the outer edge of the layer y = b(x) gives 

F dx ""^ffn |2i 4. v 1^ + i M ^ V 6u dy = 0 (2.2) 

J J V ax ay p dx by^J 

Xi 0 ^ ^ 

If Su does not depend on the choice of and Xg, 

then from the conitition of arbitrariness of Xj. and Xg 
the following ec|.uafcion will be obtained: 

h(x) 2 N 

r/Su du I dp du\ ^ ^ 

/ ( V. ^ + V -~ + - ^ "W -ZS dy = 0 ( 2.3) 

J V Sx dy p dx &y® / 

0 . 

It is seen readily that for various values of 5u 
different integral conditions will be obtained including 
the above-mentioned integral conditions of "7on Karman, 
Leibonsonj and (rolubev. Thus, for example, setting 



UAOA Tii No, 1070 



9 



6u = e = const ant 
yields, after integrating, the equation of Von Earmam: 

h h 
d r ^ . d r dp 



dx 



J pu. 4. - ^/ pu dy, = - h g - ^ Q (3.4) 



vfliere U(x) is tjie velocity at the outer edge of the 
layer. In the same manner > setting is an arliitrary, 

infinitely small quantity not depending on the coordi- 
nates) 

k 

fiU = €U 

< 

gives th.e syiBtem of integral conditions of GolulDevi 

h 

r / bvL &u 1 dp a*u\ k ' ^ /„ =N 

/ {urr-+v-- + --il-v ) uiy =0 (2.5) 

\ ox oy p dx dy*/ 

(k = 0, 1, 2, . . . .«) 

or, after some transformations: 



h . h 

d r ui-+^ • Tjl^+i d r 

/ ' ' dy.~ ■ ' ■ — / u dy 

dxj k+1 fc+ldxj 
o 0 

' / xfay - / u''-* (1)° ( 2.6) 



p dx 

0-0 

For k = 1, there is the equation of leilienson. In 
deriving the system of equations (2.6), use was made 
of the arbitrariness of the variation 6u, 

If in the system of equations (2,5) the author 
passes from the variables x and y to the variables 
X and I = u/TJ(x), which, in the region of the layer 
up to the separation point, are in a reciprocal 1: 1 
relation, the system of equations ihns o'bta;ined:' 
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L (x, i)i de = 0 (k» 0, 1, 2, «,) 

o 

may he considered as the conditions of setting equal to zero 
all the moments of the continuous function L(x, |), Thus 
if possible to find a velocity distrihution in the "boundary 
layer which satisfies the infinite system of equation (2,5) 
or (2,6) such profile would also satisfy the Prandtl equation- 
that is, would give an exact solution of the problem of the 
boundary layer, 

The question as to convergence of the method in using 
a finite number of equation (2^5) requires, of course, in- 
dependent investigation in dach c^ase. If such convergence 
exists the condition of Von Karman would be the equation of 
the f i^rst^ approximat ion, the combination of the equations of 
Von Karman and Leibonson would correspond to the second approx- 
imation, and so forth. 

The equation (2,3) may be considered as the integral 
condition corresponding to the general principle of Jourdain, 

In this equation the velocity variation plays the part of the 
"w/jigjat" of the integral mean. In the integral method of Von 
Karman the weight is equal to unity, each element of the in- 
tegral being assigned the same weight. In the Leibenson equa- 
tion the velocity u serves as the weight. Since the velocity 
u becomes zero at the wall and assumes the maximum value at 
the edge of the layer it may be expected that the velocity pro- 
file satisfying the single condition of Leibenson will approach 
nearer the true value at the outer region of the layer than 
nea^r the wall. The same is true with regard to the conditions 
of aolubev for k = 1, 

lilide use may be made of the arbitrariness in the choice 
of the variation 8u for the purpose of obtaining any partic- 
ular weight of the integral mean. 

The equations of Leibenson— CrOlubev, written in the form 
(2.5) and (2,6), as the equation of Von Karman in the form 
4'2~,-^4i)-,y m&y hsf ..a^^ a layer with finite thickness. 

In using the more accuraVe theory of t'Kl '"ilyMii^^^^^^^ 
these equations 3,ose their validity because the .integrals be- 
come infinite. ^■^pther more convenient form of the equations 
may be * indicated, valid for both the finite and the asymptotic 
layers. 
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- equations it is necessary to pass, in 

the equations of the "bouhdary layer, from the velocity' . -a 
to the so-called "velocity deficiency" q = U - u where U 
is the velocity at the outer edge assumed in advance as a 
given function of x. Use is made of the continuity equation 
to ohtain 

Multiplying t>oth sides by q^ and integrat ing with 
respect to y between the limits zero to the thickness of 
the boundary h or to infinity in both cases after simple 
transformations and integration by parts there is obtained: 



00 , h 



= (k+l)V J q^ |,§ dy (2.7) 



By analogy with the so-called "displacement" thickness' 
6* and "momentum thickness" 6** equal, respectively, to 



00 



h 



, 00, h 



are introduced the thicknesses 
00 h - 




k A . vk 



0 ' 0' 

co.h 00, h (hic** hjc*- h£+i ) 
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Equation (2,7) then assumes the form 

^c4.« ^^^^^ ^k+i ^ (k+DU^+^Tjrh*^^ 

CO . 



= (k+1) V r q^ ^ dy 

By* 



0 

Dividing hpth sides of this equation "by U^"*"^ finally gives: 



dh* * „ . 



ij- j^(k+2) hj^:, + (k+i)hj^j 



dx 

.h 



05 I"* 

2 



/(-^) i^®^'' '^-^^ 



This system of equations corresponds to the general 
equation (2.3) for 8u = (U— u) € = q^ e and, in contrast 

with the equations of Leihens on— Soluhev, may be applied "both 
to the layer of finite thickness and the asymptotic layer, 

Por k = 0 the thicknesses hf and h^** "become the 

usual displacement thickness 6* and momentum thickness 
6** a.nd equation (S,8) passes over into the equation of Von 
Karman in the form (In this paper to denote the thickness 
of the layer the letter h, and not 6, is used so as not 
to confuse this ,sym"bol with the variation sign. In equation 
(2.9) the usual notation is kept.) which is now best taken as 

(26- 4. 6.) '. (2.9) 

dx U pU® 

where Tj, is the friction stress at the wall; 

' = ^. (2.10) 

\oy/y=5 0 ■ 



Equations (2,9) and (2,10) may be considered, respec- 
tively, as the integral and differential (with respect to 
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the coordinate . y) determinations of. the frictional stresses 
for a glveii velocity distribution u(x, y)^ The approximate 
differentiation of a given function u(x, y) is a very in- 
accurate Operation for determining the friction. In general, 
the integral method of determination (2,9) is to he preferred. 
The only exception may be the region where the velocity pro- 
files rapidly change their shape and inaccuracy in deter- 
mining the derivative d6**/dx may show up. This will "be 
the case in the region near the separation of the boundary 
layer. Thus, in determining the local frictional stresses, 
it is recommended to use the integral method (2, 9) , everyvrhere 
except in the region near the separation, Tor determining 
the separation point/ however, it is necessary to obtain 
Internal solutions — that is, solutions which are nearer the 
true values at the walls than at the edge of the layer and at 
the separation point equate to zero the friction determined 
by the differential method, 

I'rom this point of view the Von Karman method, based oh 
the application of equation (2,9), in which the right side is 
replaced by (2,10), is no other than the condition that the 
local frictional stress determined by the differential and 
integral methods be the same. This condition, . in using only 
a single parameter (thickness of the boundary layer), is dif- 
ficult to satisfy even by imposing on the velocity profile a 
large number of boundary conditions because by using this con- 
dition Use is not made of the specific advantages of the dif- 
ferential and integral expressions for the friction. But by 
use of. "a method of approximate . deterainat ion the^ve^ocity pro- 
file that does not involve the use of the Von Karman equation 
it is possible to apply equation (2,9) as the integral expres- 
sion of the local frictional stress. 

These considerations apply entirely to the Iieibenson — 
G-olubev equations, sincie the Von Karman equation is included 
and. is fundamental in the system of successive approximations 
(setting the moment of zero order equal to zero). 

As usual in the practical applications of direct methods, 
only the first approximation is of significance, since to ob- 
tain the second and following approximations is a problem of 
great^^if f i^c^Ul.ty, It v;ill be shovm how, by making use of the 
general equation (2,3) and' the "integrar expression for the 
friction, it is possible to obtain ^ fjrst approximation more 
accurate than the solution of Von Karman— Pohlhaus en, 

3., It is best to start with the simplest case of the 
flow about a plate, yellowing the idea of a boundary layer 
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of finite thickness, the velocity distr ihution in a cross- 
secf3.on of the layer ie expressed in ...the form: 

u = Vs^(r^) = (y/h) (3.1) 

where U = constant is the velocity of the approaching flow 
and So^n) is a function satisfying some condition on the 
surface of the plate (r\ « 0) and at the outer edge of the 
layer Cn = i), if the thickness of the layer h(x) is the 
fundamental and only parameter, in the velocity distribution 
eterl JheJ'' ^® ^7 varying this paraml 

6U » Ugo« (ti) 8n » - Vngo' (tl) (3.2) 

h'.-. 

Simple computations give; 

0 

Substituting these expressions in the general equation (2.3^ 
yields, because of the arbitrariness of 6h: ' 

- »ieo'= (^jT godn) dn^^y^ t,^^. g^" dn= o 



On the variation 6ti is imposed tjae condition of be- 
coming equal to zero on the boundary of the interval ti = 0. 
l=li fox this by (3.2) it is required that g« ' (l) = 0 
Then integrating by parts t-l^e second integral iS the abo-^e 
equation gives . 



dx " ■' ■ u 

where 

1 



o 



NACA TM No. 1070 



15 



For the condition h( 0) = 0 there is obtained: 

h = a J V x/U (3»4) 

The two equations (3,4) and (3,3) constitute a solution 
of the prohlem of the thickness of the boundary layer, 

The local frictional stress may be determined by the 
differentiation method by the formula: 



o-r by the integral method on the basis of (2.9)? 

0 

where $ denotes the constant: 



P = y god- So) (3.7) 



In the assumed averaging the weight Su becomes zero 
at the wall and at the outside boundary of the layer. It 
should be expected here that the velocity profile obtained 
will approach the true one in the middle part of the boundary 
layer and give sufficiently accurate values for the friction 
determined integrally by equation (3,6) 

The application of the simple parabolic distribution 

So = 2^1 - Ti^ 

satisfying the condition u = 6 for y si 0, 'u' = tJ" and 
bu/by = 0 for y3ih(x) gives, as simple computations show, 
the friction formula 



- 0.331 yfipU*/x 
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differing by only 0,3 percent of the accurate formula of 
'Blasius.^^while' application of the same di strfbution in the 
Von Karman method leads to the formula; " 



Tq * 0,365 J M-pTJ^/x 
exceeding the exact solution by 10 percent 



XABLZ 1 





a 


/uoU = 






3n - 3n* + 


3 , 25 


0,92 


0. 17 


0.39 


1 ^ - i i' 

8 2 


3.6 


.42 


.25 


0.32 


3 - 3 ^ 


3.8 


'.35 


.27 


0.31 


sin HH 
2 


4.42 


,36 


•.31 


).33 


2n - TV® 


4.95 


\40 


.331 


5,36 


2n~ 2ti' + 11* 


5,62 


.36 


.331 


J. 34 


Blasius' exact 
solution 




.332 


'.332 









If go^^) subjected to a larger number of conditions 
by making use of the values of the second and third derivatives 
at the T^'all and edge of the layer, another series of values of 
thje. ,.CQ!g;ff,icie_nts in the .integral formula for the friction may 
be obtained. ' 'fn " table 1 the result sf'^^f computations, arranged 
in the order of increasing thickness of the layer (coefficient 
a,) are given. As may be seen from the second and fourth col- 
umns, X'/ith increasing thickness of the layer (coefficient a) 
the friction determined integrally at first increases to a 
value near the true one and then decreases. In the last column 
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are given the values of the friction coefficients according 
to To'hlhausen f or which the rule jus; t given, does not hold. 
It may "be easily shown that in the particular case of a 
plate the friction coefficients according to Pohlhausen are 
the geometric mean of our differential and integral coeffi- 
cients* 

4, It is not difficult to generalize the ahove device 
of varying the thickness of the layer to the case of two— 
dimensional flow about an arbitrary body. 

It is well to start with the simplest example of a 
parabolic approximation of the velocity profile in the bound- 
ary layer, Ihis, as has been shown, immediately leads to a 
very simple formula for the solution in finite form of the 
friction problem. In contrast to the plate where the velocity 
of the outside flow is the same throughout, in the case of an 
arbitrary body, the velocity U in the formula 



(4.1) 



will be a given function TJ(x) of the abscissa x, 

The value of the variation of the velocity 8u remains 
the same as in the previous section, but the expression in 
parentheses in the fundamental equation (2,3) will becomo 
somev/hat- more complicated after substituting (4,1), This 
expression will be given later in a general form, but for 
the present it is enough to show that, in the concrete case, 
(4,1), after very simple computations, eq.uation (2,3) is 
reduced to the. differential equa-tion? 



17 h» 



13 U' 



630 h- 126 U 



V 

3h3U 



= 0 



(4.2) 



With introduction, as is customary, of the function z=h®/v 
equation (4.2) assumes the form 



^ = 7.65 i: 2 ='24.7 1 
dx U U 



(4.3) 



If the above equation is rewritten in the form 
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- 2 4.7 - 7y65 U«g 

' " 'dx' ' a 

tt is seen that if U( 0) = 0 for x s 0 the point x = 0 
is a singular point of the equation. Imposing on the deriv- 
ative dz/dx the condition that it finite for x = 0 
results in the equation: 

24,7 - 7.65 Uo«Zo =0, 

for determining th« initial value of Zq for x = 0, From 
this equation it follows that the initial value of the psw- 
raaeter o:f Pohlhausen \ in the present case is equal to 



^o = ^o'^o = 3-23 (4.4) 

Thus to the equation (2,S) there is added the initial coo- 
dition 

z = 3,23/Uo ' for ' x = 0 
Equation (4,3) is linear and has the solution: 

X 

2 = J&f 24,7 ir^»«^/ U®'®^ dx (4.5) 

0 

If, as in the case of the plate, If ^ 0 for x = 0 
the initial value of z is equal to zero - that is, the 
initial thickness of the layer is equal to zero,. 

Again making use of the integral determination of the 
friction, equation (2,^) results in (8** « 2h/15 , 6* = h/s) 



In the. case of the plate (U' = 0) formula (4,6) reduces 
to the formula previously obtained for the friction at the 



i^J:k.jm\]f:^:^^ f f'^.^ d?t) (4.£) 
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plate, (The computations in the given example for the pareu- 
■bo'lic" approximation were carreid out "by Stepaniantz. The 
latter also gives a generalization of formula (4,6) for the 
case of a flow with axial symmetry ahout a hody of revolution,) 

Bq.uation (4,6) gives a sufficiently accurate determinar- 
tion of the friction in that part of the boundary layer where 
U' is near zero. As has already been pointed out, the ac- 
curacy for the plate exceeds 0,3 percent. Near the forward 
stagnation point, however, and at the separation point where 
IT« > 0 or U« < 0 the accuracy drops, the friction obtained 
is too low at the forward stagnation point and too large near 
the separation point. Since in practice in the diffuser part 
a transition occurs from the laminar to the turbulent layer, 
this defect is not of great significance, Pormula (4,6) may 
be recommended for practical pomputations in view of its 
great simplicity, Jor the function U(x) graphically given 
the integrals may be easily computed wit^ the aid of mechanical 
devices while the derivative U'(x) is determined by graphical 
differentiation. The use of the method of Von Karman— Pohlhausen 
for a parabolic velocity profile also leads to a simple fric- 
tion forniula but one of very small accuracy (for a plate the 
error, as has already been pointed out, is eq.ual to 10 percent); 
in the diffus«r region it gives a greatly exaggerated friction 
and no separation at all while formula (4,6) gives a separation 
point though somewhat retarded, 

For a comparative estimate of formula (4,6) with the 
formula- of Yon Karman— Pohlhausen n^ar the forward critics^l 
point X ss 0, U = 0 the Von Karman— Pohlhausen formula is 
given: 



In the immediate neighborhood of the point x s 0 the 

velocity of the external flow be represented in the fo,rm 

U = cx. Then, as simple computation shows, formula (4,6) 
gives : 




0.365 U°*®^ 





(4.8) 



while formula, (4,7) leads to the expression 
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(4. 8 • ) 



In the simplie case considered TI st cx the accurate 
solution may also be given (reference 10); 




— = 1.23264 / - i; (4.9) 
pU* «/cx 

It is seen that in the iinniediate neighborhood of the 
critical point "both comparison formulas give lowered solutions 
while the order of error is about 10 percent. When U(x) 
deviates fr.bm 'the straight line U * cx formula (4.6) gives 
a solution closer to the true pne than formula (4;»7), which 
gives too great a value for the friction. Present velocity- 
profiles for flight angles of attack are such that only very 
near to x = 0 is there a region for which a straight line 
U = cx approximately holds true, Even for very small values 
of X the curve rapld,ly passes Over into a region, of almost 
constant value of the velocity- IT and therefore small TJ', 
lor such a type of curve formula (4.6) should give sufficiently 
accurate values for the friction. 

In order to estimate the accuracy of formula (4.6) in 
the region beyond the maximum velocity (diffuser region), 
consider the example of the profile 

U = bo - b^ X (4.10) 

repre-senting the simplest linear drop of the velocity of 
the external flow beyond the maximum point. In equation (4,10) 
the coefficient bo gives the maximum value of the velocity 
Uq for X * 0 and b, !s U' = constant, the slope of the 
velocity curve. Substituting for x a new variable x* = 
bj^x/bo gives 

X X* 
_ . 7. 65 

IT®'®'* dx « / (l-x*)®*®^ dx» 



O 



7.65 b, 
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Suljstituting the. valiLie of this integral in formula (4.6) 
results in 



.8 .83 9 



■^^ . 1:^ Ji:::^^^ cg- 1.35(1^ x*r^'^°] 

Let the following magnitude independent of ho and h^ 
he considered: 



it 



1— X* 



1— :..x* 



The magnitude Tj^ « / ( 1— x* ) was obtained hy Howarth 

(reference 11) as a result of an accurate solution of the 
problem for the velocity profile igiven. for comparison, . in 
tahle 3, the values of this magnitude are given as computed 
hy Howarth and hy formula (4.11), 



fABliB 2 



X* 


According to Howarth 
■ 1 - X* 


By formula (4.11) 

1 - X* 


0,0000 




00 ■ 


. 0125 


2,773 


2.78 


. 025 0 


1.817 


1,82 


, 05 00 


1.064 


1.21 


.lOpO 


,.345 


.645 


, ,1200 


.000 




0,22 




. • - . 000 



The tahle shows that only for small values of x* — 
that is, only at the start of the diffuser part is there 
sufficiently good agreement between formula (4,6) and the 
accurate solution. As the point of separation is approached 
(x*.=: 0,12), as has been point out earlier, the method becomes 
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increasingly inaccurate, Tor x* = 0, 05 the approximate 
solution exceeds tlie.^accurate almost 15 percent. Setting 
the right side of (4.11) equal to zero the abscissa of the 
separation point Xg* = 22 is found, while the accurate 
value is Xg* = 0,12, The Pohlhausen formula (4,7) gives a 
greatly exaggerated friction and leads to the incorrect con- 
clusion of the absence of separation. Thus for a value as 
small as x» = 0,0125 the increased friction, according to 
Pohlhausen, constitutes about 25 percent. 

In 1939 the new formula of Wright. and Bally (reference 
12) for the frictional stress appeared; 



This formula, attractive on account of its extreme 
simplicity, was derived by Wright and Bailey on the basis 
of an experimentally confirmed assumption of the possibility 
(for small divergence and convergence) of substituting in 
equation (2,9) for 6*/6** and i5**/6.x the corresponding 
values for the plate, Formula (4,12) is accurately true for 
a plate but contains an errpr that is immediately evident 
to the eye: namely, the friction is determined by the values 
U and U' at a given point and does not depend on the pre- 
vious development of the layer^ It is readily seen that at 
the forward part it is extremely .inaccurate. Thus if, as 
before, • TJ » cx there is obtained, by equation (4»12),. 



a result more than twice the accurate value (4,9), The 
proposed formula (4,6) is almost as simple as the formula of 
Wright and Baily and at the same time is free from the above- 
mentioned defects. 

More accurate approximations will now be considered, 

5, As is known, the velocity profile satisfying the 
boundary conditions at the wall and at the outer edge of a 
layer of finite thickness must be represented in the form: 
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u = U (x) Tso + Xgi(Ti) 3 " (5.1) 

where Tj = y/h{x) and \ the known parameter of Pohlhausen 
is given by 

\ a U« h*/ V (5.2) 



The functions Sq^^) ^^'^ Sj(^) polynomials the 

degree and coefficients of which depend oh the choice of the 
particular boundary conditions. 

The general form of relation (5.1) is Icept and the 

fundamental eq^uation (2.3) is set up, where use will be made 
of the device of varying the thickness of the layer, or what 
amounts to the same thing, varying th$ parameter X, This 
will' give the following eq^vation (& dot above a letter 
denotes differentiation with respect to t\ and a dash 
differentiation with respect to x): 



6u s U (Sgo + \Sgi + gi S\) 



= U ( go + Xg^ ) 



/ 1 SJ^^ 



+ Ug^ 6\ 



= l - i T»g,^ \ - I Hgo] 8^ (5.3) 
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The expression in parenthesis in equation (2.3) may be 
given in the form: 

= u(g, -*- [w {g,^^g,)-*- Ue, V -i^u{^ - ^) (i, -»-Xft)]- 
Removing parenthesis and making certain simplifications 
yieldd: 

Substitution of the obtained expression L(u,v) and 6u 
in equation (2.3) results in (for complete arbitrariness 
of 6X ) a differential equation with respect to X: 

where 

The coefficients a^, a^, .... a^^, ^re constants de- 
pending on the form of the functions gg and g, - that 
is, on the choice of the boundary conditions. These con- 
stants are determined, respectively, as the integrals 
with respect to r\ from 0 to 1 of the following func- 
tions of i: 
_ 1 . . • 

= - 4- ri, ^ - A jf. A-i. - 1 J -/.(/I - T 
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^ — —- J- V» ^ftft— ft] f I A j I '« I ft j ~ 

— 4" {*« ~ T ft f *" • 
5. = _ A r,ir, ^ ft' - ir. I j -H - 4" I ft ~ ft I ft '''' ) 

-•- T Jft - T (ft - T '''J ( ft f ft ft f ft '''' ) • 

= (», - 4" ^.) ^ ft' - ir. J ft A j — 5- (ft - T 'ii) if. Jft 

(1 l \ 1 

ft j ft A -♦- ir, f ft Jii j Y (ft — y 1 A ) ft j ft . 

^1 =— T *i f ft «'•»-»- 4" (ft — T ""ft) I ft J ft ft" i ft ] ' 
«n = 4 (ft ~" T " ft ) ft f ft • 

0 ^ 

At request of the author, A. P. Krol has made computa- 
tions for the boundary conditions: 

«=0. for Jf = p, 

u^U. ^=0. ^=0 for » = A. 

in satisfying which the functions go'l) and gi(n) 
haye the form: 

fc(r)=-:2ii — 2«»-*-Ti«, ft(»i) = (Ti — 3r»H-3Ti» — V)/6. (5lS> 

The values of the computed coefficients are given below: 
Q.Sn4KS7, mi a< — OilDOOOSaOS. 

«is ai0234«6S. a, = OL00OO0M7. 

a,^— aO0H)6»2. •( K — <U>n7«)4>. 

a000l3l4S. aOOKOT. 
QU00000379. ait« 0000001. 

mt=- omn«m. mooooom. 

a,.— 000127321. 

The functions xtX] and <tr(^) may be tabulated iii 
advance 
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, .. Eauation (5,6) has the singular point x = 0, \ = 

the usual investigation shows IHiat X© the least posii^ive 

root of the numerator of and is eq.ual to = 6,73,, 

The second singular point will the point x = Xj,, \ = 0 
of the minimum pressure. 

As for the prohlem of the integral determination of the 
friction t at th* wall, there is first found: 

8* = (0,300000 - O,0p8333X)h = a,(X)h 

(5,10) 

6** = (0.1175 00 - 0,001058\- 0. 000ll02\^)h = b( \)h 



Then t is determined hy formula (2,9): 



If desired, the differentiation may "be carried out in 
the first term on the right and, the derivative of \ with 
respect to x from equation (5,6) substituted. 

To estimate the relative accuracy of the method, it is 
necessary to proceed as "before. It is already known that 
for' small values of TJt the integral formula (S.ll) gives 
excellent agreement with the accurate solution of Blasius 
(difference about 0,3 percent, table 1, in which are given 
the values of Ti for the plate). 

The behavior of the friction near the forv/ard critical 
point must be considered next, It is noted that U'sc, 
TJ»« = 0, (*\/'ii'^)x«:o ~ ° TJ = cx and for sufficiently 

small X it may be assumed that neglecting powers higher 
than the first — \ s \q s 6,73, Then by equation (5,11) 

there is obtained: 

-;^='^'^o C2b(Xo) + a (^ )3 yr i = 1.18 J^i (5,13) 
ptJ^ " °vcx Vex. 

# 

a result approximately 4 percent less than the accurate 
(4,10), Thus the chosen polynomial of the fourth degree 
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satisfying conditions (5,8) determines the friction near the 
critical point with almost two— and— ar-half times the accuracy 
of_t he previously considered polynomial of the second degree 
and in the remaining region (excluding the region near the 
separation point) gives the same accuracy (of the order of 
0,3 percent). The same approximation, according to Pohlhau— 
sen, gives near the critical point: 

A ^ 1,20 /^^ ' ~ (5.13) 

that is, only 3 percent less accurate, while with increasing 
distance from the critical point the accuracy of the solution 
(5,11) considerably exceeds the accuracy of the Pohlhausen 
solution, 

The method considered in the present and preceding sec- 
tions of estimating the accuracy of the solution permits the 
conclusion that the order of the error in the forward part 
does not depend on the steepness with increasing U(x); the 
absolute error of the proposed solutions will accumulate in 
the region of increasing Tj(x) and will he smaller the 
smaller the ahscissa interval where the velocity approaches 
its maximum value. 

The computations conducted by A, 1>, Krol for the diffuser 
region, showed that the friction computed by the method just 
discussed is much nearer the true value than that computed by 
the method of Pohlhausen and other methods. Only in the im- 
mediate neighborhood of the point of separation does the method 
fail. The reasons for this fact require special investigation. 

6, The method presented still requires additional im- 

provoment and simplification. It may be mentio^ned that L, G, 
Stepaniantz generalized the method of varying tHe thickness of 
the layer to the case of a flow with axial symmetry about a 
body of revolution. The results obtained show that the method 
may be successfully applied to the practical computation of 
bodios of revolution (of airship form). 

It is hoped that later on the method of choice of the 
most suitable forms of the velocity varia:tions for a single 
and, several xarying parameters will be made more accurate. 
It is apparent that there are no db j*e"^i" 6ns In princ iple to 
the application of the method not only to the laminar but 
also to the turbulent boundary layer. 

Translation by S, Reiss, 
national Advisory Committee 
for Aerona,utics. 
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